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We exploit analytic continuation to prolongate to the region of real chemical potentials the 
(pseudo)critical lines of QCD with two degenerate staggered fermions at nonzero temperature 
and quark or isospin density obtained in the region of imaginary chemical potentials. We de- 
termine the curvatures at zero chemical potential and quantify the deviation between the cases 
of finite quark and of finite isospin chemical potential. In both circumstances deviations from a 
quadratic dependence of the pseudocritical lines on the chemical potential are clearly seen. We 
try different extrapolations and, for the nonzero isospin chemical potential, confront them with 
the results of direct Monte Carlo simulations. We also find that, as for the finite quark chemical 
potential, an imaginary isospin chemical potential can strengthen the transition till turning it into 
strong first order. 
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1. Introduction 

The determination of the QCD phase diagram in the temperature-quark density plane is be- 
coming increasingly important, due to its impact in cosmology and in the physics of compact stars 
and of heavy-ion collisions. The first-principle nonperturbative approach of discretizing QCD on a 
space-time lattice and performing numerical Monte Carlo simulations is plagued, at nonzero quark 
chemical potential, by the well-known sign problem: the fermionic determinant is complex and 
the Monte Carlo sampling becomes unfeasible. Analytic continuation is amongst the possible al- 
ternatives to solve (approximately) the sign problem [1-3]. It consists in performing Monte Carlo 
simulations at imaginary chemical potential /I = i^m'- where there is no sign problem. The results 
obtained at imaginary chemical potential are then analytically prolongated (jJ.m — s> — z'jU) at real 
values of the chemical potential. There are however limitations due to ambiguity in the interpo- 
lation and nonanalyticities and periodicity [4], so that reliable estimations are expected only for 
Re(/i)/r < 1, where T is the temperature. 

In previous works [5-7] we have studied the analytical continuation of the pseudo critical 
line in the case of SU(2) with n/ = 8 staggered fermions and finite quark density and SU(3) with 
tif = 8 staggered fermions and finite isospin density. It was found that the nonlinear terms in the 
dependence of J3 C on /i 2 in general cannot be neglected and that the extrapolation to real /i may be 
wrong otherwise. We have also studied [8] SU(3) with tif = 4 staggered fermions and finite quark 
density. In this case we observed deviations in the pseudocritical line from the linear behavior in 
jj, 2 for larger absolute values of jx 2 and we saw that there are several possible extrapolations to real 
jU that are in agreement with each other up to \ijT ~ 0.6. 

In the present study [9] we consider two-flavor QCD in presence of a quark or an isospin chem- 
ical potential in the standard staggered discretization for fermion fields, whose partition function, 
in the standard staggered discretization for the fermion fields, reads 



In Section 2 we present results on the analytic continuation of the critical line, T c {jx) from imagi- 
nary to real jU in the case of a finite isospin chemical potential jUj so , where simulations are available 
for both imaginary and real /ij so and on the analytic continuation of the quark chemical potential 
}i q . In Section 3 we make a comparison between the two theories at finite \i q or jUi so , quantifying 
systematic differences for quantities like the curvature of the pseudocritical line at zero chemical 
potential. In Section 4 we study the nature of the transition as a function of the isospin chemical 
potential. 

2. Analytic continuation of the pseudocritical line 

We performed numerical simulations on a 16 3 x 4 lattice (apart from some special cases where 
we varied the spatial size to investigate the critical behavior) for bare quark mass am = 0.05 cor- 
responding to m n ~ 400 MeV. We used the Rational Hybrid Monte Carlo (RHMC) algorithm, 
properly modified for the inclusion of quark/isospin chemical potential. Typical statistics have 
been around 10k trajectories of 1 Molecular Dynamics unit for each run, growing up to 100k tra- 
jectories for 4-5 j8 values around the pseudocritical point, for each /I 2 , in order to correctly sample 
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Figure 1: (Left) Pseudocritical couplings in finite isospin SU(3) with rif = 2. Negative values of 
(Hiso/(nT)) 2 correspond to imaginary chemical potentials. (Right) Fit to real and imaginary isospin chemi- 
cal potential data using a ratio of a fourth to second-order polynomial (ratio (4,2) fit). 



the critical behavior at the transition. The pseudocritical /3(/i 2 ) has been determined as the value 
for which the susceptibility of the (real part of the) Polyakov loop exhibits a peak. In fig. [l] (left) 
the data for the pseudocritical coupling versus (jUi so / '(nT)) 2 are shown. To interpolate these values 
we can exploit the ratio of polynomials: 



ao + a 1 (At/(7rr)) 2 + e7 2 (At/(7rr)) 4 + a 3 (At/(?rr)) t 
l+a 4 (n/(nT)) 2 



(2.1) 



The fit to all data requires (see fig. [j]) at least a ratio of fourth order to second order polynomial 
(ratio (4,2) fit) and gives a # 2 /d.o.f. = 0.6. If we consider data with (}i/nT) 2 > — 0.375 2 a linear 
(in (fi/nT) 2 ) polynomial works quite well (# 2 /d.o.f. = 0.95), contrary to our previous findings 
for other theories [5-8] where nonlinear corrections are more important for imaginary values than 
for real ones. The interpolation to only imaginary (/j.i so /(nT)) data using a atio (4,2) fit gives a 
^ 2 /d.o.f. = 0.49. We have also interpolated imaginary isospin chemical potential data using the 
implicit relation between j3 c - and jU 2 



C7 2 ( j 8 c (M 2 ))| 2 -loop = a 2 (/3 c (0))| 2 -loop 



1 +Ax + Bx 2 



l+Cx 

and the following interpolating function ("physical" fit) given in terms of the physical units x 
H/(nT) and T/T c (0), 



(2.2) 



Tc(0) 



1 +Ax + Bx 2 



(2-3) 

1 +Cx 

with T = 1/ (N t a(fi)), we also get a very good # 2 /d.o.f. = 0.53 and a corresponding prolongation 
to real values that works quite well (see fig. ||] left). Another quite good interpolation of imaginary 
isospin data is attained by means of a sixth-order constrained polynomial fit, where the coefficient 
of (jU / (nT)) 2 is fixed at the value derived from a linear (in (jU / (nT)) 2 ) fit at small imaginary chem- 
ical potential data. In fig. |2] (right) the extrapolations to real isospin chemical potentials together 
with results from simulations at real values are shown. We can see that several extrapolations agree 
up to n/(nT) < 0.2. Therefore we may conclude that different interpolations that well reproduce 
imaginary data, lead to distinct extrapolations (as we have seen [8] for n-f = 4 SU(3)). 
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Figure 2: (Left) Fit to the pseudo critical couplings in finite isospin SU(3) with «/ = 2 according to the 
"physical" fit Eq. (2.3). (Right) Extrapolation to real isospin chemical potentials of the sixth-order con- 
strained, ratio fourth to second-order polynomials and "physical" fits (only the border of the 95% CL band 
have been reported). Data points (circles) are the results of Monte Carlo simulations performed directly at 
real isospin chemical potential. 



Let us move now to the nonzero quark chemical potential simulations. In this case the sign 
problem prevents us to perform simulations at real values of the quark chemical potential. In fig. || 
(left) we can see that the ratio (4,2) interpolation used in the case of isospin chemical potential is 
well suited here too, giving a % 2 /A.o.f. = 0.60. If we tried a linear fit (in ji ) we got a largely 
unsatisfactory # 2 /d.o.f. = 2.87. As for the isospin chemical potential we also tried the "physical" 
fit (Eq. (23) to the imaginary quark chemical potential data. The result, as shown in fig. || (right), is 



good also in the present case {% /d.o.f. = 0.51). Assuming that it is possible to extrapolate down to 



T = the relation T c (n)/T c (0) versus /I (Eq. (p3|), we get the following estimate for the chemical 
quark potential at zero temperature: 




Hc{T = 0) = \ -T c (0) = 3.284(65)r c (0) , (2.4) 



to be compared with \i c (T = 0) = 2.73(58)J C (0) of ref. [10] with n f = 2 Wilson fermions. In 
fig. |] the extrapolations at real values of the quark chemical potential starting from three different 
successful interpolating functions at imaginary chemical potential values are compared. The three 
analytic continuations begin to deviates at ^l/(jtT) > 0.1. However two of these extrapolations 
(in particular the ratio of polynomials and the "physical" fit) continue to be in good agreement. 
Moreover we observe that in the case of isospin chemical potential the ratio of polynomials is pre- 
ferred, but we cannot claim this is the interpolation to use for analytic continuation since systematic 
differences between finite quark density and finite isospin QCD cannot be excluded. 



3. The curvatures of the critical lines 

To obtain the curvatures of the critical line at }i = ((t/j8 c (/x 2 )/(3f^ 2 )|^ = o) for isospin and 
quark chemical potentials respectively, we tried a common fit to all data we have collected for the 
critical couplings at quark chemical potential and isospin chemical potential (in the latter case we 
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Figure 3: Fits to the pseudocritical couplings at finite quark density: ratio of a 4th- to 2nd-order polynomial 
(left) and "physical" fit according to the function (2.2) (right). The dashed vertical line indicates the boundary 
of the first Roberge- Weiss sector, (fim) / '{nT) = 1 /3. 




Figure 4: (Left) Extrapolation to real quark chemical potentials of the 6th-order constrained, ratio (4,2) of 
polynomials and "physical" ratio fits (only the borders of the 95% CL band have been reported). (Right) 
Comparison between the extrapolations to real quark and isospin chemical potential of the fits linear in 
(fJ,/(nT)) 2 . Data points (circles) are the results of Monte Carlo simulations performed directly at real 
isospin chemical potential. 



included also data at imaginary values): 



jScG^.Miso) = &(0) + a q 2 + a iso 



(3.1) 



We included in the fit as many data as to have a reasonable # 2 /d.o.f. and we obtained: a q = 
-0.3997(87), a ho = -0.3606(67), j3 c (0) = 5.32370(57) with a £ 2 /d.o.f. = 0.93. Therefore we 
can conclude that the curvatures of the critical lines respectively for isospin and quark chemical 
potentials differ up to 4 standard deviations. By expressing the curvatures in terms of dimensionless 
quantities [11] 



with 



R 



q/iso 



Ida 



ft(0) 



1 q/iso 



/AtisoV 
\kTJ 



N c 



2J3 C (0) 3 p L (Pc(0),m q 



^q/ho ) 



(3.2) 



(3.3) 
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Figure 5: (Left) Normalized plaquette distributions at the pseudocritical coupling for different spatial lattice 
sizes and jj,- l!i0 / (nT) = 0.475/. (Right) Maxima of the plaquette susceptibility as a function of the spatial 
volume for n iso /(nT) = 0.475/. 



where p L = a^ is the 2-loop lattice beta-function, we get R q = -0.515(11) and/? iso = -0.465(9) 
in agreement with Refs. [10, 12, 13]. Therefore: 

Rq _ iso = = = 0.098(26) ~ 10% . (3.4) 

Rq Oiq 

This could be the first evidence of the &{\/N^) difference between the two theories at small chem- 
ical potential [14-17]. 



4. Order of the phase transition at imaginary chemical potentials 

The phase structure at finite T and imaginary chemical potential may be important of its own 
and teach us something about the nonperturbative properties of QCD also at zero and small chemi- 
cal potential. The phase transition at the Roberge-Weiss endpoint could in principle have influence 
also far from the endpoint. For nj = 2 the Roberge-Weiss transition is first order for small and high 
quark masses and second order for intermediate quark masses [18]. In the present study we have 
am = 0.05 so that we expect a second order Roberge-Weiss phase transition (at /Iim/ {nT) = 1/3) in 
the case of the quark chemical potential. On the other hand we expect that imaginary isospin chem- 
ical potential may strengthen the transition as an imaginary quark potential does: a first-order tran- 
sition could be manifest along the pseudocritical line (even for our quark mass value am = 0.05). 
In fig.|| (left) we display the normalized plaquette distributions at the pseudocritical coupling for 
different spatial lattice sizes (L s = 12, 16,20), while in fig.|| (right) we can see the maxima of the 
plaquette susceptibility that scale linearly with the spatial volume. Therefore we conclude that for 
nf = 2 staggered fermions of mass am = 0.05 the transition is first order at Hi so /{7iT) = 0.475/ 
and there is possibly a critical point along the line at some smaller value of jJ.i S0 /{nT). Such 
non-trivial behavior resembles what happens for quark chemical potentials [18-21] and may have 
consequences on the general structure of the QCD phase diagram. 



References 

[1] M. G. Alford, A. Kapustin, and F. Wilczek, Imaginary chemical potential and finite fermion density 
on the lattice, Phys. Rev. D59 (1999) 054502, [|hep-lat/9807039||. 



6 



Two-flavor QCD at finite quark or isospin density 



L. Cosmai 



P. de Forcrand and O. Philipsen, QCD phase diagram for small densities from simulations at 



imaginary mu, Nucl. Phys. Proc. Suppl. 119 (2003) 535-537, [ |hep-lat/0209084| ]. 

M. D'Elia and M.-P. Lombardo, Finite density QCD via imaginary chemical potential, Phys. Rev. D67 



(2003) 014505, [hep-lat/0209146] 



A. Roberge and N. Weiss, Gauge theories with imaginary chemical potential and the phases of QCD, 
Nucl. Phys. B275 (1986) 734. 

P. Cea, L. Cosmai, M. D'Elia, and A. Papa, Analytic continuation from imaginary to real chemical 



potential in two-color QCD, JHEP 02 (2007) 066, [ |hep-lat/0612018| ]. 

P. Cea, L. Cosmai, M. D'Elia, and A. Papa, The critical line from imaginary to real baryonic chemical 



potentials in two-color QCD, Phys. Rev. D77 (2008)051501, [arXiv: 0712 . 3755] 



P. Cea, L. Cosmai, M. D'Elia, C. Manneschi, and A. Papa, Analytic continuation of the critical line: 



suggestions for QCD, Phys. Rev. D80 (2009) 034501, QarXiv: 0905 . 1292j |. 

P. Cea, L. Cosmai, M. D'Elia, and A. Papa, The phase diagram of QCD with four degenerate quarks, 



Phys. Rev. D81 (2010) 094502, [arXiv: 1004 .0184] 



P. Cea, L. Cosmai, M. D'Elia, A. Papa, and F. Sanfilippo, The critical line of two-flavor QCD at finite 
isospin or baryon densities from imaginary chemical potentials, Phys. Rev. D85 (2012) 094512, 



[arXiv: 1202 .5700] 



K. Nagata and A. Nakamura, Imaginary chemical potential approach for the pseudocritical line in the 
QCD phase diagram with clover-improved Wilson fe unions, Phys. Rev. D 83 (Jun, 201 1) 1 14507. 

M. D'Elia and F. Sanfilippo, Thermodynamics of two flavor QCD from imaginary chemical potentials, 



Phys.Rev. D80 (2009) 014502, QarXiv: 0904 . 1400[ ]. 

P. de Forcrand and O. Philipsen, The QCD phase diagram for three degenerate flavors and small 



baryon density, Nucl. Phys. B673 (2003) 170-186, [ |hep-lat/0307020 



J. Kogut and D. Sinclair, The Finite temperature transition for 2-flavor lattice QCD at finite isospin 



density, Phys.Rev. D70 (2004) 094501, [ |hep-lat/0407027| ]. 

D. Toublan, A Large N(c) perspective on the QCD phase diagram, Phys. Lett. B621 (2005) 145-150, 



[ |hep-th/05010 69p . 

M. Hanada and N. Yamamoto, Universality of Phases in QCD and QCD-like Theories, JHEP 1202 



(2012) 138, [arXiv: 1103. 5480] 



A. Armoni and A. Patella, QCD With A Chemical Potential, Topology, And The 't Hooft 1/N 



Expansion, Phys.Rev. D85 (2012) 125021, QarXiv: 1204 . 2405Q . 

M. Hanada, Y. Matsuo, and N. Yamamoto, Sign problem and phase quenching in finite-density QCD: 



models, holography, and lattice, arXiv: 1205 . 1030 



C. Bonati, G. Cossu, M. D'Elia, and F. Sanfilippo, The Roberge-Weiss endpoint in Nf= 2 QCD, 



Phys.Rev. D83 (2011) 054505, QarXiv: 1011 . 4515Q . 

M. D'Elia and F. Sanfilippo, The Order of the Roberge-Weiss endpoint (finite size transition) in QCD, 



Phys.Rev. D80 (2009) 111501, [arXiv:0909. 0254] 



P. de Forcrand and O. Philipsen, Constraining the QCD phase diagram by tricritical lines at 



imaginary chemical potential, Phys.Rev.Lett. 105 (2010) 152001, QarXiv: 1004 . 3144Q . 

C. Bonati, P. de Forcrand, M. D'Elia, O. Philipsen, and F. Sanfilippo, Constraints on the two-flavor 
QCD phase diagram from imaginary chemical potential, PoS LATTICE2011 (201 1) 189, 



[arXiv: 1201 .2769] 



7 



